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THE EVOLUTION OF THE PECULIAR TYPE lA SUPERNOVA SN 2005HK OVER 400 DAYS 
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ABSTRACT 

UBVRI photometry and medium resolution optical spectroscopy of peculiar Type la supernova 
SN 2005hk are presented and analysed, covering the pre-maximum phase to around 400 days after 
explosion. The supernova is found to be underluminous compared to "normal" Type la supernovae. 
The photometric and spectroscopic evolution of SN 2005hk is remarkably similar to the peculiar 
Type la event SN 2002cx. The expansion velocity of the supernova ejecta is found to be lower than 
normal Type la events. The spectra obtained > 200 days since explosion do not show the presence 
of forbidden [Fe ii], [Fe ill] and [Co ill] lines, but are dominated by narrow, permitted Fe ii, NIR 
Ca II and Na i lines with P-Cygni profiles. Thermonuclear explosion model with Chandrasckhar mass 
ejecta and a kinetic energy smaller (E'k = 0.3 x lO^^ergs) than that of canonical Type la supernovae 
is found to well explain the observed bolometric light curve. The mass of '^^Ni synthesized in this 
explosion is O.I8M0. The early spectra are successfully modeled with this less energetic model with 
some modifications of the abundance distribution. The late spectrum is explained as a combination 
of a photospheric component and a nebular component. 

Subject headings: supernovae: general — supernovae: individual (SN 2005hk) 



1. INTRODUCTION 

An impressive homogeneity in the light curves and 
peak luminosities make Type la supernovae (SNe la) 
good candidates in the determination of the extragalactic 
distance scale. Though a majority of the observed SNe la 
belong to the "normal" type (Branch et al. 1993), a num- 
ber of studies indicate significant photometric as well as 
spectroscopic differences. For example, studies of nearby 
supernovae by Li et al. (2001) indicate that 64% SNe 
la are " normal" , while 20% are of the overluminous SN 
1991T type and 16% of the underluminous SN 1991bg 
type. The peak absolute luminosities of SNe la are well 
correlated with their immediate post-maximum decline 
rate, forming a photometric sequence from the luminous 
blue events with a relatively slow decline rate to the 
faster, red, subluminous events (Hamuy et al. 1996a, b; 
Phillips et al. 1999). 

However, there are a few SNe la that are known to 
deviate from this relation, the most notable amongst 
them being SN 2002cx. This supernova was found to 
be underluminous, but had a light curve decline rate 
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Ami5(i?) = 1.29, comparable to normal SNe (Li et al. 

2003) . The early phase (< 100 days after explosion) spec- 
tra indicate line velocities lower by a factor of 2 compared 
to those of normal SNe la (Li et al. 2003; Branch et al. 

2004) . Furthermore, the late phase (~ 250 days after 
explosion) spectra are also quite dissimilar compared to 
normal SNe la and possibly consist of P-Cygni profiles 
(Jha et al. 2006). 

Interestingly, SN 2002cx is not a unique event. Jha et 
al. (2006) and Phillips et al. (2007) have shown that SN 
2005hk shows photometric and spectroscopic behaviour 
almost identical to that of SN 2002cx. Further, Jha et al. 

(2006) have shown that SNe 2003gq and 2005P also show 
spectra similar to SN 2002cx. Spectropolarimetric obser- 
vations of SN 2005hk at the early phases by Chornock 
et al. (2006) indicate low polarization levels, indicating 
that the peculiarities of SN 2002cx-likc SNe do not re- 
sult from an extreme asphericity. Based on the analyses 
of the early phase spectra of SN 2002cx, Branch et al. 
(2004) suggest that the observed lower line velocities are 
consistent with the deflagration models of explosion. Jha 
et al. (2006) report a possible detection of O i lines in 
the late phase spectra and suggest large scale mixing in 
the central region, which is also consistent with three- 
dimensional (3D) deflagration models. Phillips et al. 

(2007) find that, qualitatively the observed hght curves 
of SN 2005hk are in reasonable agreement with model 
calculations of a 3D deflagaration model that produces 
~ 0.2 Mq of s^Ni. 

Understanding the nature of this class of SNe la is 
thus quite important for the study of homogeneity and 
heterogeneity of SNe la. It may provide a caution for the 
cosmological use of SNe la and a strong constraint to the 
explosion models. 

We present in this paper the photometric and spectro- 
scopic development of SN 2005hk over ~ 400 days since 
explosion. 
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Fig. 1. — Identification chart for SN 2005hk. The stars used as 
local standards are marked as numbers 1-11. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Photometry 

Supernova SN 2005hk was observed in UBVRI bands 
with the 2m Himalayan Chandra Telescope (HCT) of 
the Indian Astronomical Observatory (lAO), Hanlc, In- 
dia, using the Himalaya Faint Object Spectrograph Cam- 
era (HFOSC), equipped with a SITe 2k x 4k CCD. The 
central 2k x 2k pixels were used for imaging which corre- 
sponds to a field of view of 10' x 10', with a scale of 0.296 
arcscc/pixcl. The monitoring of the supernova started on 

2005 November 5 (JD 2,453,680), ~ 4 days and ~ 6 days 
after its discovery by Burket & Li (2005), and Baren- 
tine, et al. (2005), and continued until 2006 Septem- 
ber 20 (JD 2,453,999). Standard fields PG0231-I-051, 
PG1047+003, PG0942-029, PG0918+027 (Landolt 1992) 
were observed on 2005 November 24 and 2005 December 
28 under photometric conditions, and were used to cali- 
brate a sequence of secondary standards in the supernova 
field. 

The data reduction was done in the standard man- 
ner, using various tasks available within IRAF. The ob- 
served data were bias subtracted, flat-field corrected and 
cosmic-ray hits removed. Aperture photometry was per- 
formed on the standard stars, using an aperture radius 
determined using the aperture growth curve, and were 
calibrated using the average color terms and photomet- 
ric zero points determined on the individual nights. The 
UBVRI magnitudes of the secondary standards in the 
supernova field, calibrated and averaged over the two 
nights are listed in Table 1. The secondary sequence 
is shown marked in Figure 1. The magnitudes of the 
supernova and the local standards were estimated using 
the profile fitting technique, with a fitting radius equal to 
the FWHM of the stellar profile. The difference between 
the aperture and profile-fitting magnitude was obtained 
using bright standards in the supernova field and this 
correction was applied to the supcaiiova magnitude. The 
calibration of the supernova magnitude to the standard 
system was done differentially with respect to the local 
standards. 

SN 2005hk was also observed on 2005 December 26, 

2006 June 30 and November 27 with the Faint Object 



Camera and Spectrograph (FOCAS; Kashikawa ct al. 
2002) mounted on the 8.2 m Subaru Telescope in the 
B and R bands. The local standard star 11 was cali- 
brated as indicated in Table 1 using the HCT/HFOSC 
photometry of the star 2 in Figure 1 , then the B and R 
magnitudes were derived by the point spread photometry 
in the DAOPHOT package of IRAF. 
The observed supernova magnitudes are listed in Table 

2. 

2.2. Spectroscopy 

Spectroscopic monitoring of the supernova using the 
HCT began on 2005 November 4 and continued until 
2006 January 12. The journal of observations is given in 
Table 3. The HCT spectra cover the wavelength ranges 
3500-7000 A and 5200-9100 A at a spectral resolution 
of 7 A. The data reduction was carried out in the 
standard manner using the tasks available within IRAF. 
The data were bias corrected, flat-fielded and the one 
dimensional spectra extracted using the optimal extrac- 
tion method. Spectra of FcAr and FcNe lamps were 
used for wavelength calibration. The instrumental re- 
sponse curves were obtained using spectrophotometric 
standards observed on the same night and the supernova 
spectra were brought to a relative flux scale. The flux 
calibrated spectra in the two regions were combined to 
a weighted mean to give the final spectrum on a rela- 
tive flux scale. The spectra were brought to an absolute 
flux scale using zero points derived by comparing the ob- 
served flux in the spectra with the flux estimated using 
the photometric magnitudes. 

Spectra during the late phases were obtained with the 
Subaru/FOCAS on 2005 December 26, 2006 June 30 and 
November 27. A 0.8 arcsec width slit was used with 
the B300 grism, which gave a wavelength coverage of 
4700-9000 A and a spectral resolution of 11 A. The ex- 
posure times were 900 sec, 1,800 sec and 3,600 sec for 
2005 Dec 26, 2006 Jun 30 and Nov 27, respectively The 
instrumental response was corrected using data of a spec- 
trophotometric standard star (GD 153 or BD-|-28d4211). 
The flux was then scaled to he consistent with the R- 
band photometry. 

3. UBVRI LIGHT CURVES 

The early phase light curves (LCs) of SN 2005hk are 
discussed by Phillips et al. (2007) and Stanishev et al. 
(2006). We present here the LCs of this supernova dur- 
ing the early as well as the late phase (ref. Table 2). 
Figure 2 shows the LCs of SN 2005hk in the U, B, V, 
R and / bands, respectively. The epoch of maximum in 
each of these bands is estimated by a cubic spline flt to 
the points around maximum. The supernova reached a 
maximum in the B band on JD 2453685.34 with magni- 
tude of 15.91±0.03. The maximum brightness in the V, R 
and / bands occured ~ +A, -1-6.5 and -1-8.5 days, relative 
to the maximum in B band, respectively. The decline 
rate parameter in the B band estimated using our data 
is Ami5(_B) = 1.68±0.05. Our estimates of the dates 
of maximum, apparent maximum magnitudes and the 
decline rate parameter are listed in Table 4. These esti- 
mates are consistent with the values reported by Phillips 
et al. (2007) and Stanishev et al. (2006). 

The LCs of SN 2005hk are compared with those of the 
normal SNe 1994D [Ami5(B) = 1.31±0.08; Richmond et 
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Fig. 2.— The UBVRI LCs of SN 2005hk. The LCs have been 
shifted by the amount indicated in the legend. 

al. 1995] and 2003du [Arm^B = 1.04 ± 0.04; Anupama 
et al. 2005], the overluminous SN 1991T [Ami5(B) = 
0.95 ±0.05; Lha et al. 1998], the subluminous SN 1991bg 
[Ami5(B) = 1.93 ±0.08; Filippenko et al. 1992; Turatto 
et al. 1996] and the pecuhar SN 2002cx [Ami5(B) = 
1.29±0.011; Li et al. 2003, Ami5(B) = 1.7±0.1; Phillips 
et al. 2007] (Figures 3-4). The early LCs (Figure 3) in- 
dicate that the pre-maximum brightening of SN 2005hk 
was faster than the luminous SN 1991T but compara- 
ble to the normal SNe 1994D and 2003du. Following 
the maximum, during the initial 20 days, the decline 
in B was faster than in normal SNe, but similar to SN 
2002cx, while the decline in V was similar to the normal 
SN 1994D. Subsequently, beyond 30 days past maximum, 
the decline in B and V was slower than all the other SNe. 
The secondary peak seen in the R and / bands of SNe la 
was found to be absent in SN 2005hk, similar to the sub- 
luminous SN 1991bg. However, beyond ^ 50 days past 
maximum the R LC of SN 2005hk had a decline similar 
to SN 1991T, which is slower than that seen in the nor- 
mal and subluminous SNe (Figure 4). The evolution of 
the U light curve of SN 2005hk was found to be faster 
than both the normal and the luminous types. 

The early phase evolution of SN 2005hk is very sim- 
ilar to that of SN 2002cx. Both objects have similar 
Ami5(B). SN 2005hk has a Ami5(B) of 1.68, while SN 
2002cx has a value of 1.7, as indicated by the the revised 
photometry by Phillips et al. (2007). On the other hand, 
at epochs beyond 15 days past maximum, SN 2005hk 
declines somewhat slower than SN 2002cx. The decline 
rate during ~ 20 — 50 days past B maximum, in the B, 
V, R and / bands are, respectively, 0.021, 0.027, 0.031 
and 0.031 mag day~^, while the corresponding decline 
rates in SN 2002cx during the same phase are 0.037, 
0.035, 0.041 and 0.035. It may be noted here that no 
galaxy template subtraction is done for the data pre- 
sented here, and the SN magnitudes may be affected by 
the host galaxy at the later phases. However, as sug- 
gested by Phillips et al.(2007), neither K correction nor 
the galaxy background contamination is sufficient to ex- 
plain the discrepancies seen in the LCs of SN 2005hk and 
SN 2002CX. 
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Fig. 3.— Early phase (< 50 days) LC of SN 2005hk together with 
those of SNe 2002cx, 1991bg, 1991T, 2003du and 1994D. The ordi- 
nate in each panel is the magnitude below the respective maximum, 
and the abscissae represent the days since respective B maximum. 
For clarity, the LCs have been shifted by the amount indicated in 
the legend. 
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Fig. 4.— Late phase (> 50 days) LC of SN 2005hk together with 
those of SNe 2002cx, 1991bg, 1991T, 2003du and 1994D. The ordi- 
nate in each panel is the magnitude below the respective maximum, 
and the abscissae represent the days since respective B maximum. 
For clarity, the LCs have been shifted by the amount indicated in 
the legend. 



4. REDDENING AND COLOR CURVES 

The reddening within our Galaxy in the direction of 
SN 2005hk is E{B - V)q^x = 0.022 (Schlegel et al. 1998). 
The reddening within the host galaxy may be estimated 
from the methods suggested by Phillips et al. (1999), 
Ahavilla et al. (2004) and Lha (1996). The Phillips et 
al. (1999) and the Altavilla et al. (2004) methods give 
E(B—V^Yiost = 0.20 and 0.09, respectively, while the Lira 
(1996) method indicates a higher value. Spectropolari- 
metric observations of SN 2005hk (Chornock et al. 2006) 
indicate an interstellar polarization of 0.27% produced 
by the host galaxy which, for the standard dust polar- 
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Fig. 5.— The U-B, B-V,V-R and V-I color curves 
of SN 2005hk. Also shown for comparison are the color curves 
of SNe 2002CX, 1991T, 1991bg, 2003du and 1994D. The abscissae 
correspond to days since respective B maximum. 

ization efRciency, corresponds to E{B — V^)Host ~ 0.09. 
Owing to the peculiar nature of SN 2005hk, the methods 
based on normal SNe la may not be applicable. Hence, 
adopting the estimate for the host galaxy extinction as 
indicated by the spectropolarimetric observations, we es- 
timate the total reddening towards SN 2GG5hk to be 
E{B - V)tot = 0.11. 

The {U - B), {B -V), {V - R) and {V - I) color 
curves of SN 2005hk are shown in Figure 5, compared 
with other SNe. The color curves of all SNe have been 
corrected for reddening using the Cardelli et al. (1989) 
extinction law and the E{B — V) values of E{B ~V) — 
0.13 for SN 1991T (Phillips et al. 1992), E{B - V) = 
0.034 for SN 2002cx (Li et al. 2003), E{B - V) = 0.04 
for SN 1994D (Richmond et al. 1995), E{B -V) ^ 0.05 
for SN 1991bg (Turatto et al. 1996) and E{B-V) ^ 0.01 
for SN 2003du (Anupama et al. 2005). 

The color curves of SN 2005hk are very similar to SN 
2002cx. The {U — B) color evolution is very similar to 
SN 1991T. The {B - V) and {V - R) color evolution 
of SN 2005hk in the pre-maximum and maximum phase 
are very similar to the other SNe la except SN 1991bg, 
which had a considerably red color. Beyond 5 days, 
until ~ 20 days after B maximum, both {B — V) and 
{V — R) colors get increasingly redder compared to other 
SNe, but are still bluer compared to SN 1991bg. The 
{V — I) color on the other hand gets redder than other 
SNe 5 days before maximum, and also does not show 
the dip seen in the color curves of other SNe la at ~ 10 
days post-maximum. Beyond ~20 days after maximum, 
the {V - R) and {V - I) colors of SN 2005hk follow the 
trend of other supernovae and become marginally blue, 
while {B — V) continues to be redder than other SNe, 
including SN 1991bg. 

Interestingly, despite a very similar color evolution dur- 
ing the pre-maximum phase to ~ 20 days after B max- 
imum, SN 2005hk appears to be marginally bluer than 
SN 2002cx at later phases. 

5. BOLOMETRIC LIGHT CURVE 



5.1. Behavior of bolometric light curve 

The bolometric LC of SN 2005hk is estimated using 
the UBVRI photometry presented here along with the 
NIR YJH photometry reported by Phillips et al. (2007). 
The magnitudes were reddening corrected using the value 
estimated in the previous section and the Cardelli et 
al. (1989) extinction law. The corrected magnitudes 
were then converted to monochromatic fluxes using the 
zero points from Bessell et al. (1998). The bolomet- 
ric fluxes were derived by fitting a spline curve to the 
U,B,V, R, I, Y, J and H fluxes and integrating over the 
wavelength range 3100A to 1.63/J,m. In the later phases 
when only B, V, Ror V,R magnitudes were available, the 
bolometric flux were estimated by applying a bolometric 
correction to the available magnitudes. Since SN 2005hk 
belongs to a rare class of SNe la, the bolometric correc- 
tion is derived from the light curve of SN 2005hk itself 
rather than assuming the corrections based on normal 
SNe la (Suntzeff 1996; Contrado et al. 2000). The bolo- 
metric corrections were estimated based on the last three 
points in the LC for which the BVRIJH bolometric flux 
is available. A simple average of the bolometric correc- 
tions for these days give values of B.C = —0.926 ±0.089 
in B, B.C = 0.294±0.045 in V and B.C = 0.677±0.027 
in R band. 

The bolometric LC is plotted in Figure 6 (red circles). 
Also plotted in the figure are the bolometric LCs of the 
normal (but slightly under luminous) la SN 1992 A (black 
squares), the subluminous la SN 1991bg (blue triangles) 
and SN 2002cx (green triangles). The peak luminosity 
(il/boi = —17.7) is fainter than the canonical value for 
normal SNe la. It suggests that a smaller amount of 
^^Ni is synthesized during the explosion. 

The decline of the bolometric LC after the maximum 
is slower than that of SN 1992A although they have a 
similar maximum luminosity, as also shown by Phillips 
et al. (2007). Our observations extend the bolometric 
LC to 400 days since explosion. At such late epochs, 
the bolometric luminosity of SN 2005hk is still brighter 
than that of SN 1992A. In addition, the difference in the 
luminosity between the two SNe becomes large, reaching 
^ 1 magnitude at > 250 days after the explosion. The 
bright late phase luminosity indicates a more efRcient 
trapping of the 7-rays from decaying '^^Co in SN 2005hk 
compared to SN 1992A. 

5.2. Light curve modeling 

Phillips et al. (2007) suggest that qualitatively the ob- 
served properties of SN 2005hk are consistent with 3D 
deflagration model (e.g., Ropke et al. 2006, Blinnikov 
et al. 2006), showing the bolometric and multi-color LC 
model until ^ 80 days after the explosion. Branch et 
al. (2004) and Jha et al. (2006) also suggest a 3D de- 
flagration model for SN 2002cx. Here we give the first 
investigation of the bolometric LC of this class of SNe 
until ^ 400 days after the explosion. 

We use an LTE radiation transfer code (Iwamoto et 
al. 2000) for the computation of the bolometric LCs. 
For the 7-ray transport, a gray atmosphere is assumed. 
For optical radiation transport, electron-scattering and 
line opacity are taken into account. The electron- 
scattering opacity is evaluated by solving the Saha equa- 
tion while the line opacity is crudely assumed to be con- 
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Fig. 6. — Bolometric LC of SN 2005hk (red circles) compared 
with a synthetic LC with the EOS model (red line). Black squares 
and blue triangles are bolometric LCs of SNe 1992A and 1991bg 
(Suntzeff 1996), respectively. Model LCs for these two SNe are 
shown in black and blue lines, respectively (based on W7 but differ- 
ent mass of ^''Ni, see Table 5). Green triangles are the bolometric 
LC of SN 2002CX. In order to fix the rise time of the LCs of SNe 
1992A and 1991bg, B band rise time is assumed as trise = 19.5 X s 
days (Conley et al. 2006), where s is stretch factor of B band LC 
(Goldhaber et al. 2001). The B band rise times of SNe 1992A 
(s = 0.80) and 1991bg (s = 0.68) are estimated as 15.6 and 13.2 
days, respectively. Since the LC of SN 2005hk is deviated from the 
normal sequence of SNe la, we simply assume trise = 17 days as 
in Phillips et al. (2007). The magenta line shows the synthetic LC 
with the E008 model. 

stant (0.1 cm^ g"-'^; Iwamoto 1997) for simphcity. We 
use the W7 deflagration model (Nomoto et al. 1984) as a 
standard density structure and abundance distributions. 

First, the mass of ^^Ni [M (^^Ni)] is varied from the W7 
value (Af(^^Ni) = 0.59Mq) to obtain a good agreement 
with the maximum luminosity of SNe la shown in Figure 
6. In this procedure, the total mass of iron group ele- 
ments is kept constant to conserve the explosion energy. 
The black and blue lines in Figure 6 are the synthetic 
LCs with A/(56Ni) = 0.17 and 0.088Afo, respectively. 
Both at the early phase (ref. inset in Fig. 6) and the late 
phase, the LCs are in good agreement with SNe 1992A 
and 1991bg (Suntzeff 1996), respectively ^. However, the 
model sequence with various Af (^^Ni) never reproduces 
the brightness of SN 2005hk at > 50 days after the ex- 
plosion. 

As also suggested by Phillips et al. (2007), we now 
try to explain the LC of SN 2005hk with less energetic 
models to fit the LC around maximum. We construct 
a less energetic model simply by scaling the structure 
of W7 homologously. The mass of the burned material 
should be adjusted accordingly by taking nuclear energy 
production into account. In our model, the mass fraction 
of each burned element is reduced by the same fraction 
[{Ek + i^B)/(£'w7 + Eb)] in every radial sheU. Here Ek, 
Eb and E'vvy are the kinetic energy of the SN ejecta, 
the binding energy of a progenitor WD (r^ 0.5 x 10^^ 
ergs), and the kinetic energy of W7 model (^ 1.3 x 10^^ 
ergs), respectively. The mass of the burned elements 
reduced by the procedure above is compensated by un- 
burned C and O, with the mass fraction of X(C)=0.5 and 

^ This shows that our crude assumption of constant line opacity 
(0.1 cm^ g~^) is reasonable. 
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7. — The density structure of model EOS compared with that 



of W7 (1 day after the explosion). The kinetic energy of model EOS 
is _Ek = 0.3 X lO^^ergs. The density structure is scaled from that of 
W7 homologously. This makes the density at the inner part higher 
than the original. 

X(O)=0.5. 

The red line in Figure 6 is the synthetic LC with 
A4^(56Ni) = O.18A^0 and Ek = 0.3 x lO^^ergs. This 
kinetic energy is similar to that of the 3D deflagration 
model used in Phillips et al. (2007). The density struc- 
ture of this model is given in Figure 7 together with that 
of the original W7. This model gives a reasonable fit to 
SN 2005hk at late phase as well as at early phase. In the 
figure, we assume B maximum of SN 2005hk occurs 17 
days after the explosion, which gives the best fit of the 
LC (see also Phiflips et al. 2007 ^°). The mass of ^^Ni de- 
rived from the model is somewhat smaller than O.24A^0 
derived by Phillips et al. (2007). It may be caused by 
a different distribution of ^^Ni in the ejecta. Since ^^Ni 
in the model in Phillips et al. (2007) is extended to the 
outer layers, a part of ^^Ni does not contribute to the op- 
tical light. In this paper, we call our scaled, less energetic 
model E03. 

The behavior of the model E03 LC is understood by 
simple analysis. If the opacity is assumed to be constant, 
the time scale of the bolometric LC {tlc) is scaled as 



(Arnett 1982), where Afoj is the mass 



)\^3/4 1^-1/4 

TLC oc M^f E^ 
of the ejecta. Therefore, if Aioj is fixed (Chandrasekhar 
mass), the lower energy gives a slow decline after the 
maximum. The late phase evolution of LC is mainly 
determined by the optical depth to 7-rays. Since the 
optical depth is higher for lower energetic ejecta (r oc 
Af2.£;-\ e.g., Maeda et al. 2003), the LC of model E03 
declines slowly at late phase. 

The characteristic parameters of the LC models are 
summarized in Table 5. In the lower energetic model, 
unburned C and O should be abundant, being ^ O.8Af0 
in total. The presence of large amount of unburned C-fO 
is also seen in 3D deflagration model (Ropke et al. 2006), 
which is used in Phillips et al. (2007). Although our 
model is very simple, it mimics the properties of the 
3D deflagration model. Our model is further verifled by 
spectral modeling in §6.3. 

^'^ The rise time of 15 it 1 days is estimated by Phillips et al. 
(2007), using Sit upper limit at 15 days before B maximum. Our 
model LC assuming the rise time of 17 days is consistent with the 
upper limit. 
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Fig. 8. — Pre-maximum spectral evolution of SN 2005hk. The 
phases marked are relative to date of B maximum. The spectra 
have been corrected for the redshift of the host galaxy, but not 
corrected for the reddening. For clarity the spectra have been 
shifted vertically. The telluric lines have not been removed from 
the spectra and are marked with the symbol ©. 

6. THE SPECTRUM 
6.1. Evolution during early phases 

The spectral evolution of SN 2005hk during the early 
phases are shown in Figure 8 - Figure 12. The pre- 
maximum spectra (ref. Fig. 8) show a blue continuum, 
dominated by Fe III lines, with weak Si ii and Ca ii H&K 
absorption (see §6.3.1). Comparing the pre-maximum 
spectrum with other normal Type la as well as peculiar 
Type la (Fig. 9), it is seen that the spectrum of SN 
2005hk is very similar to that of SN 1991T. However, 
the minima of the Ca ii H& K, Fe ii, Fe III and Si ii 
lines, indicate that SN 2005hk has much lower expansion 
velocities that decrease from ^ 6900 km s^^ on day —6 
to ~ 6200 km s^^ on day —4. The similarity of the pre- 
maximum spectrum of SN 2005hk with SN 1991T-hke 
events has also been noted by Phillips et al. (2007) & 
Stanishev et al. (2006). 

The post-maximum evolution of the spectrum of SN 
2005hk is presented in Figures 10 and 12. The spectrum 
and its evolution are quite different from normal SNe la, 
as can be seen from Figure 11. The spectrum and its 
evolution closely matches that of SN 2002cx (Li et al. 
2003, Branch et al. 2004, Jha et al. 2006). The spectrum 
is dominated by Fe il and lines due to Co li, Na i and Ca ii 
are also clearly present (see §6.3.1). Some Cr li lines may 
contribute to the absorption around 4800 A (Branch et 
al. 2004). The photospheric velocity decreases to ~ 4000 
km s~^ by day -1-10. The flux in the blue continuum is 
found to decline steadily, and by day -1-14, the continuum 
is quite weak below 4500A. This could also be an effect of 
line blanketing due to Fe ii lines. The spectrum remains 
almost unchanged from day -f 10 to day -1-24. 

The spectrum of day -1-62 is very similar to that of day 
+4:5, except for the narrowing of lines (ref. Fig. 12). The 
hues in the 6500-6800 A, identified with Co ii (Branch et 
al. 2004), which by day -f32 started weakening steadily 
with time (but they could be C ii, see §6.3.1). 
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Fig. 9. — Pre-maximum spectrum of SN 2005hk compared with 
those of Type la SNe 1991T and 2003du at similar epoch. The 
phases marked are relative to date of B maximum. The spectra 
have been corrected for the redshift of the host galaxy, but not 
corrected for the reddening. For clarity the spectra have been 
shifted vertically. The telluric lines have not been removed from 
the spectra of SNe 2005hk and 1991T, they arc marked with the 
symbol ©. 




Fig. 10. — Spectral evolution of SN 2005hk from phase +3 days 
to +14 days relative to date of B maximum. The spectra have been 
corrected for the redshift of the host galaxy, but not corrected for 
the reddening. For clarity the spectra have been shifted vertically. 
The telluric lines have not been removed from the spectra and are 
marked with the symbol ©. Line identification is based on Branch 
et al. (2004) 

6.2. The late phases 

The spectra of SN 2005hk were obtained during the 
late phase on days -1-228 and -1-377 (Fig. 13). As in 
the case of the early phase spectra, the spectrum of SN 
2005hk during these phases is very different from the 
spectra of normal SNe la at similar phases. The spec- 
trum of SN 2005hk is compared with those of other nor- 
mal and peculiar Type la supernovae in Figure 14. The 
difference between the spectrum of SN 2005hk and spec- 
tra of other Type la supernovae is obvious. While the 
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Fig. 11. — Post-maximum spectrum of SN 2005hk compared with 
those of Type la SNe 1991T and 2003du at similar epochs. The 
phases marked are relative to date of B maximum. The spectra 
have been corrected for the redshift of the host galaxy, but not 
corrected for the reddening. For clarity the spectra have been 
shifted vertically. The telluric lines have not been removed from 
these spectra and are marked with the symbol ©. 
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Fig. 12. — Spectral evolution of SN 2005hk from phase +21 days 
to +62 days relative to date of B maximum. The spectra have been 
corrected for the redshift of the host galaxy, but not corrected for 
the reddening. For clarity the spectra have been shifted vertically. 
The telluric lines have not been removed from the spectra of SN 
2005hk taken with HCT and are marked with the symbol ©. The 
line identification is based on Branch et al. (2004). 

other SNe la show blends of strong forbidden lines due to 
iron [Fe ii] (5159A), [Fe ill] (4500-4800 A) and forbidden 
lines due to cobalt [Co iii] (5890A) (Kuchner et al. 1994) 
in their late phase spectrum, no signature of these lines 
is seen in the spectrum of SN 2005hk. 

The late time spectrum of SN 2005hk is dominated by 
Fe II lines. The Ca li NIR lines are strong and Na I is 
also present. The spectra of days -1-228 and -1-377 (ref. 
Fig. 13) also show forbidden lines due to [Ca il] 7291, 
7234 A and [Fe ii] at 7155 A and 7453 A. The ratio of 
the fluxes of Ca ii NIR and the [Ca ii] lines on day -f-228 
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Fig. 13. — Late phase spectrum of SN 2005hk at +228 days and 
+377 days relative to date of B maximum. The spectra have been 
corrected for the redshift of the host galaxy, but not corrected for 
the reddening. For clarity the spectra have been shifted vertically. 
The telluric lines have been removed from the spectrum. The line 
identification is based on Jha et al. (2006). 

indicate densities of ~ 10^ cm^"^ in the line emitting re- 
gion. By day +377, the [Ca il] lines are stronger than the 
Ca II NIR lines, and the line ratios imply temperatures 
< 4500 K and density ^ 10* cm^"^ (Fransson & Chevalier 
1989). The mass ejected during the explosion can also be 
estimated using the number density arrived at with the 
Ca line ratios. Assuming that the ejecta is moving with 
a constant velocity of 1000 km s~^ over a period of ^ 
380 days after the explosion and a mean atomic weight 
A = 40, the estimated mass of the ejecta is ~ O.25M0. 
As pointed out by Jha et al. (2006) the mass of the ejecta 
arrived at is indicative as it does not take into account 
factors such as clumping or more complicated structures 
in the line emitting regions. It may however be noted 
that a high concentration of the mass in the central part 
is also suggested by the spectral modeling (§6.3). 

The shift in the central wavelength of the [Ca ii] for- 
bidden emission lines indicates an expansion velocity of 

500 km s~^ on day +228 and ~ 300 km s~^ on day 
+377. The fuU width half maximum (FWHM) of the 
emission lines indicate a velocity dispersion of ^ 1000 
km s~^ on day -1-228 and ^ 500 km s^^ on day -f377. 
Similarly, the emission component of Ca ii NIR 8542 A 
line indicates an expansion velocity of 235 km s^^ on 
day -f 228 and ~ 210 km s~^ on day -1-377, with velocity 
dispersions ^ 2100 km s^^ and ~ 1200 km s~^, respec- 
tively. Low velocity O i at 7002 A and 7773 A absorption 
features are clearly identified in the late phase spectrum. 
These features were identified, for the first time in any 
SNe la, by Jha et al. (2006) in the "nebular" spectrum 
of SN 2002cx. The velocity of these lines, measured at 
the minimum is 800 km s-i. Though there IS an over- 
all similarity between the spectra of SN 2002cx and SN 
2005hk at late phases, there are some noteworthy dif- 
ferences also. Comparing the spectrum of SN 2005hk on 
day +228 with the +227 day spectrum of SN 2002cx (Jha 
et al. 2006), it is found that the expansion velocities and 
the velocity dispersions are higher in SN 2005hk. The ve- 
locity dispersion based on the [Ca ii] and the Ca ii lines 
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Fig. 14. — Late phase spectrum of SN 2005hk at +228 days and 
+377 days relative to date of B maximum compared with those 
of other Type la SNe 2003du, 2002er, 1998bu, 1991T and 1991bg. 
The spectra have been corrected for the redshift of the host galaxy, 
but not corrected for the reddening. For clarity the spectra have 
been shifted vertically. 

are in fact almost twice that estimated in SN 2002cx. 

6.3. Spectral modeling 
6.3.1. Early Phases 

The observed spectra of SN 2005hk are further inves- 
tigated by spectral modeling. For this purpose, we com- 
pute synthetic spectra using the Monte Carlo spectrum 
synthesis code developed by Mazzali & Lucy (1993), 
Lucy (1999) and Mazzah (2000). We create model spec- 
tra at 5 epochs and compare them with the observed 
spectra (—6, -1-3, -1-14, +24 and -1-45 days from B maxi- 
mum). The rise time is assumed to be 17 days, as in the 
LC modeling (§5.2). 

The Monte Carlo spectrum synthesis code assumes a 
spherical, sharply defined photosphere in the ejecta as an 
inner boundary. By tracing rays of a number of photon 
packets emitted from the photosphere, the temperature 
structure in the SN atmosphere is computed through the 
flux at each radial point (see eq. (1) and (2) in Mazzali 
& Lucy 1993). Radiative equilibrium is assumed for the 
temperature determination. For the interaction between 
photons and the SN ejecta, electron scattering and line 
scattering are taken into account. Sobolev approxima- 
tion is used for line scattering. Sobolev optical depth is 
evaluated by computing ionization fractions and electron 
populations of each ion for a given temperature structure. 
Details of the physical assumptions adopted in the code 
are summarized in Mazzali & Lucy (1993), Lucy (1999) 
and Mazzah (2000). 

The input parameters of the code are the position (i.e., 
velocity, thanks to the homologous expansion) of the pho- 
tosphere (wph) and the bolometric luminosity (L). The 
parameter L is the emergent luminosity. Since some pho- 
tons emitted from the photosphere are back scattered 
into the photosphere, the luminosity at the photosphere 
(total photon luminosity emitted from the photosphere), 
is higher than the input bolometric luminosity L (Maz- 
zali & Lucy 1993). As a density structure we use the less 
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Fig. 15. — Comparison between the observed early phase spec- 
tra of SN 2005hk at (-6, +3, +14, +24 and +45 days from B 
maximum, from top to bottom) and the synthetic spectra. Lines 
that make a major contribution to the absorptions are marked. 

energetic model that gives the best agreement with the 
observed LC (model E03, see §5 and Table 5). 

The mass fractions of elements are also required for the 
computation of the synthetic spectra. We treat the mass 
fractions as parameters rather than using the abundance 
distribution in the original LC model. The mass frac- 
tions of elements are assumed to be homogeneous above 
the photosphere. We determine the mass fraction of each 
element to give a good fit of the observed line depths. 
As time goes by, the photosphere recedes and the ab- 
sorption lines are formed in the inner part of the ejecta. 
Nevertheless, we always assume a homogeneous abun- 
dance distribution at each epoch, although we change 
the mass fractions of the elements. The mass fractions 
derived by the modeling of each spectrum are thought 
to be a representative abundance ratio near the photo- 
sphere because the absorption lines are formed near the 
photosphere most efficiently. 

The parameters in the best fitting cases are summa- 
rized in Table 6. Figure 15 shows the comparison be- 
tween observed spectra and synthetic spectra. Synthetic 
spectra are reddened with E{B — V)Gai — 0.022 and 
E{B - V")Host = 0.09 (§4.). The distance modulus 
/i — 33.46 is assumed. In the figure, the flux and models 
at epoch are scaled and shifted for clarity. 

The spectrum of SN 2005hk at day —6 is in good agree- 
ment with a synthetic spectrum computed with logL 
(erg s~^) = 42.56 and Vph = 6500 km s"-'^. The strong 
absorption features around 4300A and 5000A are mainly 
due to Fe ill (Mazzali et al. 1995). The photospheric 
temperature is computed as T = 17000 K by taking into 
account the back scattering effect (Mazzali & Lucy 1993). 
This high temperature makes the Fe ill lines strong and 
also reproduces the color of the observed spectrum. 

The temperature at the photosphere is relatively 
higher than that of normal Type la SNe at similar epoch, 
and it is comparable to the temperature at very early 
phase (~ —10 days; Tanaka et al. 2007). The weakness 
of Si II suggests that the outermost ejecta of SN 2005hk 
do not contain much Si (X(Si) ^ 0.02) as in normal SNe 
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la (X(Si) > 0.5). The S i hnes are similarly weak. The 
mass fraction of S is estimated as X(S) ~ 0.02. The 
Ca II lines are also very weak and there is no clear evi- 
dence of high velocity absorptions that is usually seen at 
prc-maximum epochs (Mazzali et al. 2005; Tanaka ct al. 
2006). No C lines are seen, which gives an upper limit of 
the C mass fraction X(C) < 0.001. The rest of the ejecta 
consist of oxygen, whose mass fraction is estimated as 
X(0) ~ 0.86. 

The day +3 spectrum matches a synthetic spectrum 
with log(L) erg 8'^= 42.66 and Vph = 6000 km s^^. The 
temperature at the photosphere decreases (T ^ 9, OOOK) 
from pre-maximum epochs. The low temperature allows 
some singly ionized ions such as Si II, Ca II and Fe ii to 
appear in the spectrum. 

The feature at 6300 A mainly consists of the Si ii 
A6355, and it becomes stronger from pre-maximum 
epochs (Figs. 8 and 10). After the maximum, Fe lines 
begin to contribute to this feature (Figs. 10 and 12). No 
C lines are visible, giving an upper limit of C mass frac- 
tion X(C) < 0.001. The mass fractions of Si, Ca, Fe 
and Ni can be estimated by the visible lines (Table 6). 
The mass fraction of Si, stable Fe and ^^Ni required for 
the fitting of this spectrum are larger than those for pre- 
maximum spectra. Note again that we assume spatially 
homogeneous abundance distribution for each epoch. 

After maximum brightness, the evolution of the spec- 
tral shape is not significant (Fig. 10). Almost all the ab- 
sorption features are identified and they are mostly O i, 
Ca II, Fe II and Co ii lines. The spectra are reproduced 
by almost the same element fractions (Table 6). Estimat- 
ing the mass fraction of stable Fe separately with ^®Ni 
becomes difficult because more than 25% of ^^Ni, which 
is one of the dominant elements, has decayed into ^^Fe 
at > 40 days after the explosion. 

An interesting feature is a weak absorption at 6400 A 
in the spectra at day -|-14 (and possibly in day -|-24). 
Branch et al. (2004) suggested that this is a Co ii line. 
However, if we reproduce this feature by Co, the line 
blocking at < 4000 A is too strong, making the flux there 
too low. This line could be reproduced by the C ii A6578 
line. This identification is, however, quite marginal. If 
wc assume this feature is the C II, it requires X(C) ~ 
0.01. This is larger than the upper limit derived in the 
preceding spectra, which traces the abundances in the 
outer layer. 

The O I A7774 line is always visible in the spectra from 

day -|-14 to -|- 45. The O i line is very weak in the ear- 
lier spectra, as a result of a severe suppression of O i 
fraction due to the high temperature. In fact, the mod- 
eling suggests a high O mass fraction in the outermost 
layers (X(0) > 0.75). These facts clearly indicate that 
unburned O exists in every place of the ejecta at v > 1000 
km s~^. The mass of O derived from spectral modeling 
is discussed in §6.3.2. 

6.3.2. Late phase 

The late phase spectra of SN 2005hk consist of a com- 
bination of P-Cygni profiles of permitted lines and emis- 
sion features of forbidden lines (Fig. 13). We note that 
no strong [O I] line is seen in the observed spectra. The 
presence of P-Cygni features suggests that the ejecta are 
not completely transparent even at such late epochs. 



We study the late phase spectra using the model EOS, 
which has similar properties to those of 3D deflagration 
model. Considering the observed spectra, there seems to 
be two major problems in this model. 

(1) The ejecta of the E03 model become optically thin 
at the late phases. The photospheric velocity becomes 
zero at 250 days after the explosion if a constant line 
opacity (0.1 cm^g^^) is assumed. 

(2) In the E03 model, the mass of oxygen is as much as 
0.44 Mq. Such a model would give conspicuous emission 
lines of [O i] (e.g., Kozma et al. 2005). 

First, we create synthetic spectra with the Monte Carlo 
code used in the modeling of early phase spectra. Here 
we assume the presence of the region emitting continuum 
light at the velocity of Vph = 250 km s~^. This could be 
realized if the central part of the ejecta is denser than 
that of the E03 model. A synthetic spectrum computed 
with log (L) (erg s~^) = 40.86 and the element mass 
fraction same as at day -1-45 (blue line in Fig. 16) is in 
reasonable agreement with the observed spectrum at day 
-1-228 except for emission features at 7300 A ([Ca ii]) and 
8600 A (Ca ii). The temperature at the photosphere is 
as low as 4400 K. Most of the absorption lines are Fe ii 
but some lines of neutral atoms such as Na i and Fe i are 
also found. 

An absorption of the O i A7774 line is seen as discussed 
in §6.2 (also reported in SN 2002cx, Jha et al. 2006). In 
the synthetic spectrum, we could get the O i line with 
the mass fraction X(0) ^ 0.5. Therefore, we conclude 
that the unburned O is present in whole ejecta, down to 
very low velocities, suggesting a nearly completely mixed 
abundance distribution. This property is consistent with 
3D deflagration models. 

In modeling the early and late phase spectra, the mass 
fractions of elements are optimized assuming homoge- 
neous abundance distribution above the photosphere at 
each epoch. If we assiimc these mass fractions are char- 
acteristic values at each velocity range of the ejecta, the 
mass of each element can be computed by integrating the 
density proflle of the E03 model. In Table 5, the masses 
of C, O, Si and Fe group elements derived from spectral 
modeling are shown. Compared with the masses of the 
E03 model, the smaller amount of C and Si, and larger 
amount of O are preferred by the spectra. The nuclear 
energy release corresponding to the element abundances 
obtained from the spectral fitting is ~ 0.8 x 10^^ erg. 
Since the binding energy of the WD is ^ 0.5 x 10'''^erg, 
this leads to a kinetic energy of ^ 0.3 x 10"''^ erg, which 
is consistent with the E03 model. 

In the computation above, the cooling via lines is not 
considered, i.e., no emission lines would appear. In order 
to investigate the emission lines, we compute the late 
phase spectra using a nebular synthesis code and con- 
sider the combination with the synthetic photospheric 
spectrum (as in Mazzali et al. 2004). 

The nebular code computes the deposition from 7-rays 
with a gray opacity (0.025 cm^ g""'^) by Monte Carlo 
method. Positrons are assumed to be trapped in situ. To 
take into account the different channels in which the de- 
posited energy is consumed, we assume ineb = /neb-^7,e+ 
where L^^e+ is the luminosity deposited by the 7-rays and 
positrons from O.I8M0 of '"^^Co (originally '"^^Ni). The 
fraction, /neb) of is assumed to be balanced by cool- 
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Fig. 16.— Spectrum of SN 2005hk at +228 days with respect 
to date of B maximum (black line, 3 X 10"^'^ erg cm~^ s~^ 
added to the flux for presentation), as compared to the synthetic 
"combined" spectrum (red line, the same amount added to the 
flux) . The synthetic spectrum is a combination of the photospheric 
(blue line) and the nebular (green line) components. 

ing via nebular emission lines, while the other fraction 
is by photospheric, allowed transitions as modeled by 
the early phase spectrum synthesis code. In other word, 
the luminosity (l-/neb) ^7,e+ is deposited below the pho- 
tosphere, while the remaining luminosity (/neb)i7,e+ is 
deposited into the optically thin layer above the pho- 
tosphere. With this deposition luminosity (with vary- 
ing /neb)) non-LTE rate equations are solved assuming 
the ionization-recombination equilibrium (Mazzali et al. 
2001). Note that the code assumes that the ejecta are 
optically thin. For the calculation, we assume homoge- 
neous distribution of elements and ^^Co. The masses 
of these elements are set to be consistent with the EOS 
model, i.e., Mi^'^Co) = O.ISMq and M(0) = 0.4Mq. 

Figure 16 shows the combination of the photospheric 
(as described in the first part of this section, with the 
luminosity corresponding to M{^^Co) = 0.27Mo ^^) and 
nebular components with /neb = 0.5. We find that the 
[O i]AA6300, 6363 doublet is not strong in the model 
spectra. Indeed, the resulting [O i] is so weak that it 
does not contradict with the non-detection of the [O i], 
if /neb ^0-5, as is consistent with the observed, allowed 
line-dominated spectrum. 

This is an outcome of (a) the high density of the ejecta 
and of (b) the assumption that the elements, especially O 
and Fe, are homogeneously, thus microscopically, mixed 
together with one another. With these two conditions, 
the deposited luminosity is effectively consumed by co- 
pious [Fe li] lines rather than [O i] The high density 
is a direct consequence of the low energy SN la model. 
However, the complete mixing is a rather extreme as- 
sumption, so that we need more detailed modeling of the 
late phase spectra of SN 2005hk. 

With the nebular component, we obtain better fits 

This luminosity is higher than that of the EOS model because 
the late spectrum taken with Subaru is calibrated against the pho- 
tometry that is brighter than the model LC (Fig. 6) at 250 days 
after the explosion. 

Indeed, [Fe ll] around 5200A is too strong in the model spec- 
trum, thus smaller value of /neb is favored. 



compared to purely photospheric component, in the 
following aspects: (i) Strong emission features be- 
low - 5500A, (ii) [Fe ii]A7155 and [Fe ii]A7172, (iii) 
[Ca ii]AA7291, 7324, and (iv) Ca ii IR and [C i]A8727. 

The synthetic nebular spectrum for the E03 model re- 
sults in broader lines than the observed nebular emis- 
sions. Typical example is [Ca ii]AA7291,7324. The 
two components are separated in the observations, but 
blended in the model. This indicates that the nebular 
line emitting region is indeed more centrally concentrated 
than the model. This may again suggest the central part 
being denser than the E03 model. 

In summary, the major problems (1) and (2) above 
could be solved if the ejecta are denser than in the orig- 
inal model and the elements are microscopically mixed. 
Then, the ejecta at late phases consist of three parts: 
the optically thick innermost part {v < 250 km s^^), the 
optically thin region emitting forbidden lines {v < 700 

km s~^, i.e., the width of the [Ca ii] lines), and the outer 
layer which is thin even to 7-rays. 

7. DISCUSSION 

We have studied the bolometric LC and early and late 
phase spectra by computing synthetic LC and spectra. 
We construct a simple, less energetic model by scaling 
the kinetic energy of W7 down to 0.3 x 10^^ erg. This 
model (E03) explains the LC of SN 2005hk from early 
to late phase. The ejecta of the model contain as much 
as ^ O.8M0 of unburned C-fO material. In this sense, 
our less energetic model is similar to the 3D deflagration 
model. 

Branch et al. (2004), Jha et al. (2006) and PhiUips et 
al. (2007) suggested that the properties of SN 2002cx- 
like objects are similar to those of the 3D deflagration 
model. We confirm this suggestion by modeling the LC 
until late phase for the first time. We also find that 
the early phase spectra of SN 2005hk are also explained 
by lower energetic model. The O I absorption is always 
seen until late phase, suggesting that the elements in 
the ejecta are mixed almost homogeneously. This is also 
consistent with the 3D deflagration model. 

The late phase spectra of SN 2005hk do not show 
strong [O i] line and do show P-cygni profiles of per- 
mitted lines. They seem inconsistent with what 3D de- 
flagration model predicts, at a first look. To study this in 
detail, we modeled the late phase spectra of SN 2005hk 
by a combination of the photospheric and nebular com- 
ponents. We find that the innermost region of the ejecta 
should be denser and more centrally concentrated than 
that of our model in order to make an optically thick re- 
gion and a smaller emitting region of nebular lines. The 
absence of the [O i] line is explained if the elements are 
microscopically mixed. 

In summary, our less energetic explosion model or 3D 
deflagration model seems to be compatible with the late 
phase spectrum of SN 2005hk (and of SN 2002cx-hke 
objects) (1) if the density of the central region is higher 
than that of our scaled model and (2) if the elements 
are microscopically mixed, although it is a question how 
such extreme conditions are realized. 

The nature of this class of objects can be studied fur- 
ther if observations become available at phases later than 
presented here. SN 2005hk is still not completely trans- 
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parent at 400 days after the explosion. Only when the 
ejecta become completely thin, useful information on the 
structure of the innermost ejecta can be obtained. 

If there is a sub-class of SNe la with a lower kinetic 
energy than in normal SNe la, it could be harmful for 
cosmological use of SNe la. To show the effect, we con- 
struct a model less energetic than model EOS by scaling 
the energy of W7 model to 0.08 x 10^^ erg 1/4 of 
the kinetic energy of the EOS model). If the deflagra- 
tion is very weak and slow, such a situation will be re- 
alized. Since the nuclear energy release is only slightly 
higher than the binding energy of the WD, such an ex- 
plosion could be followed after a pulsation (pulsating de- 
flagration model, e.g., Nomoto, Sugimoto & Neo 1976). 
Since the mass of ^^Ni synthesized by such an explosion 
is thought to be smaller than that of the EOS model, we 
set Af(56Ni) = O.OQAfg (1/2 of EOS). We call this model 
E008, and the parameters are summarized in Table 5. 

In Figure 6, the synthetic LC with the E008 model is 
shown as a magenta line. The LC has a fainter peak 
magnitude than that of the EOS model because of the 
smaller amount of '"'^Ni. However, the decline rate of 
the LC after maximum is slower than in the EOS model 
because of the lower kinetic energy of the ejecta. This 
trend is opposite to the well-known LC-width relation in 
SNe la, i.e., brighter SNe la declines more slowly. 

Such explosions can be distinguished by the late phase 
LC. Since less energy leads to an efficient trapping of 
7-rays at late phases, the decline rate of the LC at > 
100 days after the explosion is slower. This is clearly 
seen in the synthetic LC (Figure 6). Although the peak 
magnitude of the E008 model is fainter than that of the 
EOS by ~ 1 mag, the late phase luminosities of these two 
models are comparable. 

A recently discovered supernova, SN 2007qd may be- 
long to this sub-class (SN 2002cx-like; Bassett et al. 
2007, Goobar et al. 2007). This supernova was reported 
to have reached a fainter maximum than SN 2002cx and 
SN 2005hk and lower line velocity, being roughly half 
that of 2002cx or 2005hk. Such properties are realized if 
the kinetic energy of the ejecta is about quarter that of 
SN 2002cx (similar to the E008 model), because of the 
relation of EkOCv"^. If this is true, the LC of SN 2007qd 
will evolve very slowly as shown in the magenta line in 
Figure 6. 

8. CONCLUSIONS 

Photometric and spectrosopic data on the peculiar 
SN 2005hk are presented. The B band light curve of 
SN 2005hk shows comparable pre-maximum brighten- 
ing, faster decline in the initial ~20 days past maximum 
and slower decline beyond ^^50 days past maximum, as 
compared to normal Type la supernovae. The fainter 
peak bolometric luminosity indicates synthesis of small 
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amount of ^®Ni in the explosion. Low expansion veloc- 
ity of the ejecta together with fainter peak luminosity is 
explained by an explosion with lower kinetic energy. A 
reasonable fit to the bolometric light curve of SN 2005hk 
is achieved with a less energetic (0.30 x lO^^erg) model 
which synthesized O.I8M0 of ^^Ni. The light curve evo- 
lution is similar to SN 20C)2cx. 

The pre-maximum spectrum of SN 2005hk is similar to 
that of SN 1991T, with much lower expansion velocity. 
The spectral evolution of SN 2005hk is very similar to 
SN 2002CX. The spectrum of SN 2005hk at late phases 
(> 200 days) is similar to SN 2002cx, except for higher 
expansion velocities and higher velocity dispersions. The 
presence of P-Cygni profiles in the late phase spectrum 
of SN 2005hk indicates that the ejecta have not become 
optically thin till our last observation. Modeling of the 
pre-maximum spectra of SN 2005hk indicates a relatively 
higher temperature which makes Fe ill lines strong. The 
presence of weak O I line at A7774 at almost all epochs 
is modeled as a consequence of high abundance of com- 
pletely mixed unburned oxygen in the ejecta. The late 
phase spectra of SN 2005hk are modeled as a combina- 
tion of the photospheric and nebular components, with 
the nebular line emitting region being more centrally con- 
centrated than is expected in the case of a lower energetic 
model. 
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TABLE 1 

Magnitudes for the sequence of secondary standard stars in the field of SN 2005hk. 



ID 



U 



V 



R 



I 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



15.896 
17.257 
19.175 

16.824 
17.842 
17.094 
17.834 
19.134 
16.366 
19.516 



± 0.020 
± 0.025 
± 0.035 
± 0.025 
± 0.027 
± 0.022 
± 0.032 
± 0.040 
± 0.024 
± 0.042 



15.425 
16.859 
17.696 

16.954 
17.645 
16.870 
17.213 
18.194 
16.470 
18.162 
18.059 



± 0.006 
± 0.007 
± 0.002 

±0.010 
±0.010 
±0.016 
±0.011 
± 0.020 
±0.007 
± 0.020 
± 0.008 



14.643 
16.067 
16.375 

16.449 
16.902 
16.150 
16.210 
17.175 
16.115 
16.944 



± 0.004 
± 0.005 
± 0.009 

±0.016 
± 0.009 
±0.010 
±0.012 
±0.014 
±0.007 
±0.010 



14.174 
15.579 
15.564 
16.097 
16.440 
15.706 
15.646 
16.571 
15.871 
16.201 
16.782 



±0.012 
±0.010 
± 0.007 

±0.016 
±0.018 
±0.017 
± 0.022 
±0.013 
±0.012 
±0.021 
±0.016 



13.777 
15.149 
14.886 
15.750 
16.003 
15.275 
15.108 
16.055 
15.623 
15.553 



± 0.005 
±0.010 
±0.010 
± 0.009 
± 0.008 
± 0.005 
± 0.004 
± 0.004 
±0.002 
±0.001 
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TABLE 2 

Photometric observations of SN 2005hk 



Date 


J.D. 
2453000+ 


Phase 
(days) 


U 






B 








V 








R 










I 








680 06 
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16 


1 OR 







014 
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_i_ 





018 


16 


UZO 


_L 

ni 





014 


16 


060 


-|- 





032 


06/ n /2f)05 


681.04 


-4.25 




_l_ 

HI 


0.046 


16 


UOD 






012 


16 




_|_ 





oil 


15 


Q'iA 
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021 


15 


979 


-t 





016 


07 /n /2no5 


682 03 


-3.27 


15.181 


-l- 
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16 
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± 
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08/1 1 /2no5 


683 06 


-2.25 


1 ^1 ^ 
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_L 
Hz 
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-L 
ZL 
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yuu 


It 
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7S/I 


It 
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15 
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± 





024 
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1 ^ 741 
1 . (41 


-X- 
31 
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16 
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_L 
3; 
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4- 
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15 
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15 
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± 
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ZC 
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U_ 
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16 222 
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1 1 y 


-L 





020 


15 


■J UO 
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16 
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-|- 
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15 


840 
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15 
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Oil 


15 
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17 
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16 
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15 


649 







015 


15 
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-L 
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18 
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17 
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16 
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± 
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16 
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009 
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31.32 








18 
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17 
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16 
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16 
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-t 
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17/12/2005 
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18 
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it 
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17 
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± 
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± 
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17 
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17 
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16 
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17 
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18 
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17 
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17 
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16 
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28/12/2005 


733.07 


47.17 
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014 


17 
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013 


17 
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it 





008 


16 
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it 





008 
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748.06 


61.99 








19 


204 


it 





031 


17 


920 


it 





014 


17 


501 


it 





Oil 


17 


142 


it 





027 


14/01/2006 


750.05 


63.96 








19 


325 


it 





045 


18 


038 


it 





035 


17 


548 


it 





027 


17 


203 


it 





015 


30/06/2006 


917.04 


230.03 








21 


578 


it 





033 












20 


140 


it 





022 












04/07/2006 


921.40 


233.31 


















20 


940 


it 





044 






















18/07/2006 


935.40 


247.14 




























20 


474 


± 





036 












20/09/2006 


999.39 


310.38 


















22 


733 


it 





168 


21 


633 


± 





099 












27/11/2006 


1066.73 


376.98 








24 


206 


it 





124 












21 


776 


it 





036 













Observed phase with respect to the epoch of maximum in B band (JD 2453685.34), corrected for the (1-tz) time-dilation factor using 
the host galaxy redshift z = 0.0118 



TABLE 3 

Log of spectroscopic observations of SN 2005hk 



Date 


J.D. 


Phase* 


Range 


Telescope 




2453000-t 


(days) 


A 





04/11/05 


679.10 


-6.17 


3500-7000; 5200-9100 


HCT 


05/11/05 


680.06 


-5.22 


3500-7000; 5200-9100 


HOT 


06/11/05 


681.06 


-4.23 


3500-7000; 5200-9100 


HCT 


13/11/05 


688.06 


2.69 


3500-7000; 5200-9100 


HCT 


15/11/05 


690.25 


4.85 


3500-7000; 5200-9100 


HCT 


17/11/05 


692.10 


6.68 


3500-7000; 5200-9100 


HCT 


20/11/05 


695.27 


9.91 


3500-7000; 5200-9100 


HCT 


23/11/05 


698.19 


12.70 


3500-7000; 5200-9100 


HCT 


24/11/05 


699.18 


13.68 


3500-7000; 5200-9100 


HCT 


01/12/05 


706.09 


20.50 


3500-7000; 5200-9100 


HCT 


05/12/05 


710.08 


24.45 


3500-7000; 5200-9100 


HCT 


13/12/05 


718.02 


32.30 


3500-7000; 5200-9100 


HCT 


26/12/05 


730.81 


44.94 


4000-9000 


Subaru 


12/01/06 


748.09 


62.02 


3500-7000; 5200-9100 


HCT 


30/06/06 


916.5 


228.46 


4700-9000 


Subaru 


27/11/06 


1067.05 


377.26 


4700-9000 


Subaru 



Observed phase with respect to the epoch of maximum in B band (JD 2453685.34), corrected for the (IH-z) time-dilation factor using 
the host galaxy redshift z = 0.0118 
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TABLE 4 

Photometric parameters for SN 2005hk 



Data 



V 



R 



epoch of max* 
magnitude at max 
Ami6(-B) 

colors at B max** 



685.34 ± 0.4 
15.91 ± 0.03 
1.68 ± 0.05 



689.49 ± 0.8 
15.71 ±0.04 



B-V 
-0.03 ± 0.04 



691.78 ± 0.2 
15.55 ±0.02 



V -R 
0.08 ± 0.03 



693.93 ± 0.3 
15.43 ± 0.03 



R-I 
-0.02 ± 0.06 



decline rate days 20-45 (mag day ^) 
decline rate days 230-380 (mag day~^) 



0.021 



0.027 
0.015 



0.031 
0.011 



0.031 



*JD 2453000+ 
colors are corrected for reddening E{B — V)total = O-H 



TABLE 5 

Characteristic Quantities of the Models 



Model 


Ek " 


M(56Ni) b 


M(C) " 


M(0) 


M(Si) ° 


M(Fe group) ^ 


W7 (92A) 


1.2 


0.17 


0.038 


0.12 


0.16 


0.95 


W7 (91bg) 


1.2 


0.088 


0.038 


0.12 


0.16 


0.95 


E03 


0.30 


0.18 


0.36 


0.44 


0.074 


0.45 


E008 


0.08 


0.09 


0.44 


0.53 


0.053 


0.32 


Spectra 






< 0.011 


0.84 


0.023 


0.44 



'^Kinetic energy of the explosion (10®^ erg s ^) 

''Total mass of ''''Ni (A/q) 

'^Mass of C contained in the ejecta (Mq) 

''Mass of O contained in the ejecta (Mq) 

"Mass of Si contained in the ejecta (Mq) 

^Total mass of Fe group elements including ^^Ni (Mq) 



TABLE 6 
Parameters of Spectral fitting 



Epoch 


fph 




log(L) d 


X(C) 


X(0) 


X(Si) 


X(S) 


X(Ca) 


X(Fe) ° 


A:(56Ni)f 


-6 


6500 


1.11 


42.56 


< 0.001 


0.86 


0.002 


0.021 


0.001 


0.020 


0.035 


+3 


6000 


1.05 


42.66 


< 0.001 


0.77 


0.020 


0.002 


0.001 


0.020 


0.14 


+14 


3500 


0.583 


42.44 


0.010 (?) 


0.55 


0.020 


0.002 


0.004 


0.010 


0.39 


+24 


1500 


0.128 


42.19 


0.010 (?) 


0.53 


0.020 


0.022 


0.004 




0.40 s 


+45 


1000 


4.90 xlO-^ 


41.88 




0.53 


0.020 


0.022 


0.004 




0.40 s 


+228 


250 


1.18x10-3 


040.86 




0.53 


0.020 


0.022 


0.004 




0.40 s 



^Days from B maximum 

''Photospheric velocity (km s" '^) 

'^Mass coordinate at the photosphere (Mq) 

''Bolometric luminosity (erg s" '^) in logarithm 

°Mass fraction of stable Fe 

'Mass fraction of radioactive ^^Ni 

SThe sum of stable Fe and ^^Ni 



